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i.  nnaoocnm  j 

!  Aa  a  part  of  tba  of  fort  to  find  a 

suitable  raproaoBtatioo  of  topography  la  a  spec¬ 
tral  model  of  tba  global  atmosphere,  a  study  baa 
bean  oada  of  tba  iapaeta  of  quadrature  aai  trun¬ 
cation  la  tba  traaaforaatiooa  between  tba  phyaloal 
iand  apaetral  deoalaa.  dpoeifleally,  rbeaboldal 
truncations  of  dlfforoat  mm  noabora  aad  two 
quadratures,  trapaaoldal  aad  Gauss-Xagendre 
(6-1),  boot  boon  Oboaaa  aa  tba  aaln  oh  jacta  of 


the  baalc  oouroa  of  laforoatloa  on  tba 
aurfaoa  topography  eonalat a  of  a  oat  of  balpht 
values  on  tba  2.5*-lntarval  latituda-longituda 
ooordinataa  furnlahad  fey  tba  U.S.  National 
Meteorological  Cantor  (NMC)  aa  a  part  of  tba 
fixed  flald  data  la  tba  KG1  mol  ZXX-A  data  aat. 

Three  measures  of  dlf  farancaa  have  boon 
oaad  to  cbaraetarlao  various  aapaota  of  tba 
Iapaeta  of  quadrature  and  truncation.  T*o, 
designated  fey  l1  aad  By  are  global  root-mean- 
aquaro  dif farancaa  daflaad  la  tba  physical  dooaln. 
If  la  the  conventional  error  of  syatheais  and 


iNsasuras  tba  eoablnad  lapaet  of  quadrature  aad 
truncation.  By.  an  tba  other  hand,  naaauraa 
tba  error  incurred  la  eeqpletlag  a  full  cycle  of 
tba  transformation  at  a  fixed  truncation  aad 
represents  tba  error  eolely  due  to  quadrature. 

qhe  third  aaasure,  I),  la  an  equivalent 
of  By  la  the  apaetral  dooaln  and  will  fee  ident¬ 
ical  la  the  total  magnitude  with  By  whan  tba 
apaetral  transform  la  exactly  invertible,  the 
ooapoeition  of  By  reveals  tba  spectral  distrib¬ 
ution  of  tba  error  of  transformation  due  to 
quadrature. 

9»  definitions  of  these  measures  and 
the  procedures  encountered  In  a  full  cycle  of 
transformation  arm  illustrated  in  Figures  1,  a 
and  B.  Additionally,  In  order  to  separata  the 
effect  of  tba  required  interpolation  from  the 
original  data  to  the  Gaussian  latitudes  from 
that  of  the  6-L  quadrature,  similar  measures  of 
differences  are  defined  with  reference  to  the 
estimates  on  the  Gaussian  latitudes  and  are 
denoted  by  priaas  la  Figure  ). 


NaamoMa 


c»(n,  to 


OHTSCAL 


ffe(M4.7» 


FIG-  1.  TRAPEZOIDAL  QUADRATURE 
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Global  circulation  model. **  Spectral  representation 

Topography 


IS.  ABSTRACT  (Continue  on  ittertt  if  neeemery  and  identify  by  Mock  lunttr)  '  —■ 

^Various  ways  of  representing  topography  that  is  a  part  of  a  global  spectral  model 
have  been  examined  using  the  notions  of  error-in-synthesis  and  error-in-reproduction 
in  measuring  accuracy  and  fidelity. 

The  trapezoidal  and  Gauss -Legend  re  quadratures  are  employed  in  conjunction  with 
different  truncations  in  the  mmerlcal  computation  of  the  transformations. 

The  results  have  led  to  the  following  conclusions: 

(1)  The  model  terrain  is  best  defined  by  the  set  of,  spectral  coefficients  in 
a  truncation  range  that  is  come  ns  urate  with  the  model  resolution.* 

(2)  The  model  terrain  in  the  real  domain  should  be  defined  on  the  Gaussian 
latltules  as  the  transform  of  the  spectral  coefficients  usinp'the  Gauss-Legendre- 
quadraturej 

(3)  Smoothing  by  a  spectral  truncation  alone  is  not  sufficient  to  insure  the 

accuracy  either  in  synthesis  or  in  reproduction.  '  .->->>,  >  -■>  r/,  f.->  1 
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(  At  tteory  of  linear  tiiufomtioo 

shown  that  007  oot  of  real-valued  data 

jl*Uj,  ♦  (*>.  J-1....^i  »-1 . .  My  be 

repreoaatad  through  tkt  fonrior  traoof  on 
by  oa  analytic  fonetiao 


•  m  4  . 

X<1,#>  -  I  ^!aiedle“* 


booh  that  XU,,  4.1  •  XU,.  ♦.)  at  all  j  and  1.  I 
kora,  a  o'J/2  iwhaa'J  la  dvea)  or  (J-D/2  (whan  , 
'J  is  odd)  ood  (W(oio^)  id  d  folynoaUl  of  dogroa 


J  id  odd)  odd  Qt<oio+)  id  d  folfvoaUl  of  doff 
|<X.  -  1>  la  ala#  for  oat*  ■  ouch  that  <4  aat  Jh*  | 
faro  o'  ooopleu  eon  jugate  pair.  j 


|  Aa  analytic  function  1(1,#) ,  on  th 

other  hand;  oay  ho  represented  In  taroa  of 
spherical  haraoalee, 

M  ■ 

*(!.♦)  -  I  I  ,XJ  fjjainfja1"* 

M-#|a-|o|  " 


in  which 


*5- jJ  fl^yJ^cyHy 


ij  *J(y>»!|/<y>dy  - 


for  all  o  and  y  ■  a in#  .  the  integral  in  (XI 
will  ha  referred  to  aa  the  Legendre  traaafoca. 
The  integrand  in  (3)  ia  seen  to  ho 


«4iCy)»S(y» 


a  polynoadal  in  y  of  degree 
(Lea*  1)  when  a  in  neon 

(1  *  y2)'/2  tinea  a  polyaooial 
in  y  of  degree  (L  ♦  n  *  2) 
when  o  ia  add. 


In  practice,  the  wave  aaplltndes  at  a  given 

latitude  ((  ere  firat  obtained  froo 

(XUj,  (|),  J-1 , . . . ,Ji  naing  the  foorier  tranafarn 

i  -inf, 

(LU.)  -  I  1(1,,  d.le  2  (4) 

J-1  - 


and  then  integral  (3)  la  roplaood  by  a  quadrature 
Mich  oay  ha  written  ao 


<  -  ItmWKWt 

*-1 


chare  yg  »  aiaft  and  g  ia  the  weight  aeaoeiatod 
with  yj .  The  reconstructed  topography  hs cones 


*(!.♦>  -  f"  f°  X*  pJleiaOe1*2 
o*-N>  n-|o| 


in  which  kg  and  lg  define  the  range  of  trunca¬ 
tion.  The  error  of  spectral  aynthesla,  x  X, 
ia  aeon  to  arioe  fren  two  poeeible  oonroea  —  ’ 
ia  the  quadrature  (S)  and  the  ether  ia  the 


The  pr possess  of  transforoatian  on 
illustrated  la  Figures  1,  2  and  3  have  hnaa 
applied  to  two  fielda  of  topography.  The 
oo-onlled  onsooothed  terrain  is  the  one  given 
in  the  roes  IXI-1  data  act.  The  ao-called 
anooth  terrain  haa  been  obtained  by  subjecting 
the  uneaoothad  terrain  to  a  nine-point  aonother 
twice.  The  anoother  la  a  product  of  two  3-point 
anon  them,  one  along  the  aonal  direction  and  the 
other  ia  the  nerldional  direction.  St  nay  he 
represented  by  a  linear  operator  lt<X  •  H 
defined  by  VU,m.,i,J)  «  wUil)v(a.J)  where 


w(i«i|  - 


1/2  if  i  •  1 
1/4  if  t  •  i  ♦  1 
•  otherwise- 


Such  a  anoother  haa  a  progressively  strong 
dsnping  affect  toward  short  waves  and  ooopletely 
eliminates  the  two-grid  interval  waves.  This 
ia  clearly  evident  in  Table  1  which  presents  the 
amounts  of  variance  contributed  by  various  spec¬ 
tral  ranges  in  both  the  unoaoothed  and  smoothed 
terrains.  These  spectra  have  been  obtained 
using  the  O-L  quadrature . 


TMLK  1  -  Mounts  of  variance  la  various  spectral  ranges  in  the 
oooothed  terrains  (unit,  n2) 
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_ _ fables  2-t  tM  2-1  tabulate  the  global 

teotum  squares  {IMS)  of  both  the  arror  in 
ijnthtilt  <t|)  and  tha  arror  in  reproduction  (Mj > 
(or  both  tarrain  fields  at  various  truncations 
,using  althar  tha  trapasoldal  or  C-L  quadrature 
vitti  tha  It  Gaussian  latitudas-  Tha  lb  Gaussian 
jlatltudas  fora  tha  aaallast  aat  required  for  tha 
•spectral  aodal  with  tha  rhoabeidal  30  truncation, 
itha  dlffaranoe  of  eharactarlatiea  between  If 
and  >2  is  readily  saan  in  tha  epposita  trends 
[of  tha  variation  of  magnitude  with  truncation 
{range.  If  dacraasaa  with  vidanlng  of  truncation  ! 

tango  as  sore  of  tha  apaetral  coupon  ants  la  tlia 
rlqlnal  fields  are  included.  I],  on  the  other 
jhand,  Increases  with  widening  of  -truncation  I 

: range  because  of  in crease  in  the  nuaber  of 
;  polynomials  whose  degrees  exceed  tha  highest 
:degree  resolvable  by  the  truncation. 

Both  tables  support  tits  preference  of 
the  c-L  over  trapasoldal  quadrature,  although 
;the  trapasoldal  quadrature  produces  a  slightly 
smaller  IMS  arror  in  If  than  tha  O-l  quadrature 
up  to  the  rhoaboidal  SO  truncation,  this  a  nail 
ndge  is  wore  than  coapensatad  by  diaadvantagaa 
found  in  other  aspects,  the  arror  in  the  global 
■can  increases  steadily  with  lncraaaiag  resolu¬ 
tion  whan  using  tha  trapasoldal  quadrature,  in 
contrast  to  the  relative  constancy  found  with 
the  G-L  quadrature  in  both  If  and  Xj  for  both 
tarrain  fields.  Nora  significantly,  the  IMS 
arror  in  reproduction  increases  wore  rapidly 


with  tha  trapasoldal  quadratore  beyond  tha 
rhoaboidal  40  truncation  in  the  anaaoothad 
tarrain  and  across  the  satire  reage  in  the 
aaoothed  tarrain. 

•  m 

l 

Further  support  for  favoring  tha  G-L 
quadrature  is  provided  by  table  3  which  pre¬ 
sents  if'  and  >2*  *°r  the  two  tarrain 
fields,  these  tables  compile  the  errors  of 
traneforaation,  had  tha  tarrain  fields  baas 
available  or  defined  on  tha  Gaussian'  latitudes, 
the  differences  between  the  corresponding  quantl 
bias  In  tables  2-1,  2-B  and  3  represent  the 
affect  due  to  tlm  antra  step  of  Interpolation 
fron  tha  2-S*-lntarval  latitudes  to  the  Gaussian 
latitudes  required  in  obtaining  tha  ■assures  Bf  . 
Sad  I2  using  tha  G-L  quadrature.  Zn  tens  of  BN 
this  aaounts  to  approximately  30-36  ■  in  tha 
un  smoothed  terrain  and  3  -  5  ■  In  tike  aaoothed 
tarrain  far  tlia  error  in  synthesis,  the  absence 
of  tha  interpolation  step  la  tha  calculation -of 
Bf'  brings  forth  the  collate  lavertibillty 
of  the  G-L  quadrature  aa  long  aa  tha  truncation 
range  doea  not  exceed  that  specified  by  the 
nuaber  of  Gaussian  latitudes  employed.  The  76 
Gaussian  latitudes  should  reproduce  exactly  up 
to  the  rhoaboidal  37  truncation,  beyond  which 
tha  arror  in  reproduction  should  increase  with 
further  widening  of  truncation  range.  Table  3 
bears  witness  to  these  theoretical  inferences. 

Sn  fact,  tha  values  of  I2'  at  rhoaboidal  60 
and  70  truncations  exoeed  those  of  12- 


miLE  2-1  •  Boot-aaan- squares  of  tlia  arror  in  synthesis  (if)  and  of  tha  arrow 
In  reproduction  (I2)  with  different  rhoaboidal  truncations  of  the  unsaootbed 


Error  Quadrature 

Xhoaboidal  Truncation 

1 

5  3 

4  30  40  SO 

Trapasoldal 

272.4 

G-L 

272.3 

Trapasoldal 

7.2 

Nhoaboidal  Truncation 


Quadrature 


Trapasoldal 

123.2 

G-L 

123.1 

Trapasoldal 

7.1 

;  Comparison*  of  tha  ms  In  Tanias  2-t, 

Mi  and  3  between  tha  unsmoothed  and  saoethad 
terrain  flalds  show  tha  large  contribution  made 
by  the  smoothing  operation  in  reducing  both  S3 
snd  I] • 


I  Tables  4-1  and  4-B  sgaurisa  tbs 

statistics  of  the  errors  In  reproduction  In  tba 
spectral  domain,  1}  and  (*3' )  for  both 
tarraln  fields,  Dpon  comparing  with  tba  corres-  ; 
ponding  quantities  in  tba  physical  domain,  ws  j 
JEind  the  C-l  quadrature  produces  smaller  differ-  : 
enema  between  the  two  domains  than  does  the 
trapezoidal  quadrature,  while  both  quadratures 
exhibit  similar  characteristics  in  the  variations 
of  magnitude  with  truncation  range  as  observed 
in  S}  and  (Sj'l*  j 

j  The  inf luenoe  of  the  differences  in  the 

terrain  fields  on  model  performance  was  assessed 
in  terns  of  the  global  root-naan-square  errors 


of  the  height  forecasts  of  the  mandatory  pressure 
levels*  For  this  purpose,  a  comparison  wea  made 
among  six  72  boor  forecasts  —  three  forecasts 
with  each  terrain  field  (un smoothed  and  smoothed) 
beginning  from  OOZ'on  15,  16,  and  17  'January  1S7S 
The  results  are  sunrlzed  in  Figures  4  and  5,  • 


•  Figure  4  shows  the  vertical  profile  of 
the  group  means  of  tba  global  root-msan-squares 
of  (1)  the  processing  error  at  tbs  initial  time 
and  (2)  the  forecast  errors  at  days  one,  too, 
and  three.  It  is  quit*  obvious  that  therm  is 
little  difference  in  both  the  magnitudes  and 
shapes  of  tha  profiles  as  a  result  of  differences 
la  the  terrain  fields.  It  is  also  clear  that 
the  processing  error  constitutes  a  email  fraction 
of  the  forecast  nrrors.  Mo  discemable  differ¬ 
ence  exists  in  the  forecast  errors  that  could  be 
ascribed  to  the  difference  in  processing  error* 
at  tha  initial  time. 


nsil  3  •  boot -mean-squares  of  the  error  in  synthesis  (*)' )  and  of  the  error  in 
reproduction  (!j*l  with  different  rhomboidal  truncation  in  reference  to  tba 
ven  values  on  the  Gaussian  latitudes  (unit  t  m 


unsmooth 


tnanooth 


•*  -  .2x10"®  is  considered  to  be  the  round-off  error. 


TUtt  4-A  -  The  square  roots  of  power  of  the  error  in  reproduction  in  the  spectral 
domain  (tjor  Ij')  in  the  un smoothed  terrain  (unit  i  m) 


Rtomboidal  Truncation 


Quadrature  Srror 


13 

24 

30 

40 

SO 

7.0 

12,4 

17.0 

22,4 

42.5 

Mhoaboldal  Truncation 


Trapesoldal  *3 


Brror 

13 

24 

30 

40  > 

* 

SO 

«) 

7.0 

12.0 

15. S 

22.2 

27.9 

To  probe  further  into  the  relationship 
between  the  initial  processing  errors  and  the 
forecast  errors,  the  group  Mans  and  standard 
deviations  of  the  difference!  between  the  global 
root-naan-square  errors  for  forecasts  with  the 
two  terrain  fields  were  calculated.  The  results 
are  shown  in  Figure  9  where  a  dot  represents  the 
group  Man  and  the  width  of  the  line  across  the 
dot  represents  twice  the  group  standard  devia¬ 
tion.  gars,  a  positive  value  indicates  a 
Mailer  error  with  the  saoothed  topography  and 
tries  versa,  Froa  Figure  5a  it  is  apparent  that 
that  use  of  the  saoothed  terrain  reduced  the 
processing  error  at  all  levels  except  the  top 
two  (i,e.,  50  nb  and  75  ab),  however,  no 
significant  trace  of  this  isproveaent  appeared 
in  the  forecasts.  The  differences  in  the 
forecast  errors  (Figures  5  b-d)  were  saaller 
than  the  differences  in  the  processing  errors 
and  were  such  saaller  than  the  forecast  errors 
theaselves,  Furtheraore,  these  differences 
were  not  statistically  significant. 


«.  COMCLOSXON 

On  the  basis  of  these  findings  we  have 
eoncluded  that  in  using  a  spectral  aodel  for 
aiaulating  and  predicting  the  global  circula¬ 
tion,  (1)  the  terrain  is  best  defined  on  the 
Gaussian  latitudes  of  the  forecast  aodel,  (2)  • 
the  Gausa-Iagandre  quadrature  is  better  than 
the  trapesoidal  quadrature  in  the  eoaputatien 
of  the  transforas  and,  (3)  the  saoothed  terrain 
la  preferable  as  the  aodel  terrain.  He  have 
consequently  defined  the  terrain  to  be  used 
in  high-resolution  aodel  as  the  oat  of  the 
spherical  haraonlc  coefficients  obtained  froa 
the  original  FOGE  data  by  first  passing  then 
through  the  9-point  saoothar  tales,  linearly 
interpolating  the  results  to  the  Gaussian 
latitudes,  and  then  tr ana forming  thee  into 
spectral  coefficients  at  the  rhoaboidal  30 
truncation.  The  aodel  terrain  is  thereby 
uniquely  defined  in  both  the  physical  and 
spectral  doaalns. 


Fig,  d.  Vertical  profiles  of  the  groi 
of  the  processing  arrora  at  day  0  and 
(a)  nnanoothed  topography  and  (b)  aaoi 


Beans  of  the  global  root-aean -squares 
mi  forecast  errors  at  days  1,  3,  and  Si 
lad  topography 


S  ao 


Fig.  4, a 


Fig,  d.b 


Fig.  5.  Vertical  profiles 
of  tho  die farancee  between 
diffarant  terrain  fields t 
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